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The study particle azimuthal anisotropy in heavy-ion collisions provides 
insight on the collective hydrodynamic expansion of the system and on its 
equation of state. The measurement of the elliptic flow, V2, of D mesons 
compared to that of light hadrons is expected to be sensitive to the degree 
of thermalization of charm quarks within the quark-gluon plasma. The first 
measurement of D° meson V2 with the ALICE detector at the LHC will be 
presented. The preliminary results obtained with the first Pb-Pb run at 
LHC show a hint of non-zero «2 in 2 < p t < 5 GeV/c. 

PACS numbers: 25.75.Nq 



1. Introduction 

The ALICE experiment pQ at the LHC is a dedicated heavy ion experi- 
ment with the goal of studying the quark-gluon plasma (QGP). 

In heavy-ion collisions with non-zero impact parameter the fireball ex- 
hibits an azimuthal asymmetry with respect to the reaction plane, defined 
by the azimuth of the impact parameter and the beam direction. Due to 
collective effects, this initial geometrical anisotropy evolves with time into 
an asymmetry in momentum space, due to different pressure gradients in- 
plane and out-of-plane. The final particle azimuthal distribution, therefore, 
is influenced by the initial geometry. It can be evaluated via the 2nd coef- 
ficient of the Fourier expansion of the azimuthal distribution, called elliptic 
flow (^2) [2]. The measurement of the heavy flavour particle t>2 is sensi- 
tive in particular to the degree of thermalization of the system, probing 
the interaction of heavy quarks, produced in the early stages of the nuclear 
collision, with the QGP. The heavy-flavour hadronization mechanism is also 
investigated [31 H] . 



* Presented at Strangeness in Quark Matter 2011, Cracow 
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Some models [3j predict charm elliptic flow to be smaller than that of 
light quarks at pt up to ~2 GeV/c and comparable at higher p t . So far, 
the heavy-quark flow has been measured via non-photonic electrons by the 
PHENIX experiment |10| . ALICE can contribute with a direct measurement 
through the fully reconstructed hadronic decays of D mesons, thanks to its 
vertex reconstruction, tracking and particle identification capabilities. 

The measurement of the elliptic flow is complementary to the measure- 
ment of the nuclear modification factor (Raa)- ALICE measured [11] a 
suppression in the D meson yield in Pb-Pb collisions with respect to pp, 
rescaled by the number of binary collisions, which probes a strong interac- 
tion of charm quarks with the QGP via energy loss mechanism. Although 
the systematics uncertainties are large, the D mesons suppression is at the 
same level as that of light charged hadrons. 

2. The first measurement of D° t>2 

This analysis is based on about 17 x 10 6 minimum bias events collected 
during the LHC 2010 Pb-Pb run. The trigger was defined as the logical OR 
of the VZERO, which consist of two scintillator arrays located at the two 
sides of the interaction point and on the Silicon Pixel Detector (SPD), which 
is the innermost detector in the central rapidity region. The centrality is 
determined via a Glauber- model fit of the VZERO signal amplitude [6] . 
The centrality class used for this measurement is between 30 and 50% of 
the Pb-Pb nuclear cross section. 

The first D° V2 measurement was performed with event plane methods 
and 2-particle cumulants [9], due to limited statistics. 

The analysis strategy consists in the selection of the D° — > K - 7r + candi- 
dates (described in ^2. 1 j) followed by the i>2 measurement with the methods 
mentioned above ( ^2.2|[2TH|) . The estimation of the systematic uncertainties 
is described in £ |2.4l and finally in ^3] the results are shown. 

2.1. Signal extraction 

A typical signature of the D° — > K~7r + decay is the presence of two 
tracks with opposite charges from a secondary vertex pointing the primary 
vertex and with an invariant mass around the D° meson mass (1.865 GeV). 
The necessary spatial resolution for the determination of the primary and 
secondary vertices is assured by the Inner Tracking System (ITS) and the ex- 
cellent track resolution is provided by the Time Projection Chamber (TPC) 
and the ITS. For each pair of tracks a secondary vertex is defined as their 
point of closest approach. The significance of the signal is optimized by se- 
lecting track pairs (candidates) that pass specific cuts on the decay length, 
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the angle between the D meson flight line and its reconstructed momen- 
tum, and the decay track impact parameter [7]. The identification of the 
kaon rejects additional background, mostly at low p t . The identification 
relies on the specific energy deposit dE/dx in the TPC and on the time of 
flight of the particle measured by the Time Of Flight detector (TOF) . 

2.2. Event plane methods 

The event plane methods correlate each candidate's azimuthal angle 
(4>) with the reaction plane angle rp). The reaction plane is estimated 
through the event plane *$> n determined by the so-called Q-vector which is a 
weighted sum of the azimuthal angles of all the TPC tracks in \n\ < 0.8 with 
quality requirements including 70 clusters in the TPC and a x 2 /n.d.f < 4. 
To remove autocorrelation, the tracks coming from the D° candidate under 
study are subtracted from the Q-vector. The ^ is used for this analysis: 

r. _ (T.f=oWi cos 2^ 
\T,i=o Wi sm20j 

In our case, the W{ are ^-weights determined on a run-by-run basis. The pre- 
cision on the measured event plane is limited by the finite number of recon- 
structed tracks. The resolution on ^ was estimated with the two sub-events 
method [T2] and the resulting correction factor c<j, 2 ((cos 2(^2 — ^rp))), to 
be applied to the observed (v2 = /cr^ 2 ), is 0.91. 

The three methods used to extract the D° V2 in 3 bins of pt are described 
in the following sub-sections. 

2.2.1. D° yield in two bins of azimuthal angle 

The first method consists in comparing the signal yield in bins of A(f> = 
(j> — ^2- The signal is extracted with an invariant mass analysis. Two 
ranges of A(p are defined: [— tt/4, vr/4) 1J [3tt/4, 57r/4) , which gives the in- 
plane yield (iVjjv) and [tt/A, 37r/4) 1J [5tt/4, 7tt/4) , which gives the out-of- 
plane yield (Nout)- The yield is extracted from the subtraction of the 
background below the exponential fit function from the total fit, in a region 
of 3a around the mean of the Gaussian fit (Fig. [U left panel) . Integrating 
the Fourier expansion in the two Acj) ranges, V2 bs can be expressed as: 

obs _ 71 N IN - N OUT 

V<2 ~ 4 N IN + Nout 

2.2.2. Method of the side band subtraction 

The cos 2A(ft distribution of D° candidates in three mass ranges is con- 
sidered: two side bands, which are then averaged, and the peak region, as 
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Fig. 1. Left panel: Invariant mass spectrum with shaded bands defining the side- 
bands (red and blue) and the peak region (black) and fit functions. Right panel: 
Distribution of cos 2A0 in different mass regions. The green area is the average of 
the left- and right-hand side-bands. p t bin 2 < p t < 5 GeV/c. 



shown in shaded areas in Fig. [H left panel. The distribution of the averaged 
side bands, shown by the green filled graph in Fig. [TJ right panel, normal- 
ized to the background fitting function under the D° peak, is subtracted 
from the cos 2A<^> in the peak region (solid black line) and finally the is 
computed as (cos 2A(/>) . 



2.2.3. Fit of v-2 vs mass 

Using the information of the fit to the invariant mass spectrum (Gaussian 
plus exponential), the fraction of signal (S) and background (B) over the 
total is estimated as a function of mass. The dependence of = (cos 2A</>) 
measured in small bins of mass (M) is drawn in Fig. [2j left panel, for 2 < 
Pt < 5 GeV/c. The fitting function is: 

*M ^( M) + - ,f M > (3) 



S{M) + B{M) 1 K ' S(M)+B(M) 

where v^ 9 (M) = pqM + p\ and v^ 9 = P2 are obtained from the fit free 
parameters po, pi, p2- 



2.3. Cumulants 

A 2-particle cumulants, the Q-cumulants method [8] was used. It was 
not possible to use 4-particle cumulants due to limited statistics. The 2- 
particle cumulants are calculated as QC{2} = (e 21 ^ 1- ^ 2 ) \ and are related 
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Fig. 2. Left panel: vf s vs Ktt invariant mass fitted with Eq. G2 ( 42.2.3|) . Right 
panel: Elliptic flow measured with 2-particle cumulants as a function of mass. The 
red closed circles correspond to the side bands mass bins and the blue open circle 
to the peak region 



to the elliptic flow by V2{2} = QC{2}. In Fig. [2j right panel, showing 
V2{2} vs mass for 2 < p t < 5 GeV/c, the closed markers are interpolated 
with a line and the background flow extracted (v^ 9 ) is subtracted from the 
measured t>2 in the mass peak region (v^ a ^) as follows: 

sgn _ S{M) + B{M) peak B(M) bkg 
V2 ~ S(M) V2 ~^{M) V2 W 

where S, B are respectively the signal and background obtained from the fit 
to the invariant mass distribution in 3a around the measured D° mass. 



2.4- Systematic uncertainties 

A 2.5% systematic uncertainty due to the different event plane resolu- 
tion within the centrality class is accounted for all event plane methods. 
For the method of the yield extraction in two A<p bins a ~ 20% uncertainty 
is induced by the signal extraction procedure. A systematic ranging from 
20% up to a factor 2-3 is estimated from the cut variation. The method was 
validated using charged tracks and comparing to official ALICE results. 
The systematic uncertainty for the method that uses side band subtraction, 
estimated by varying the width of the side bands and their binning ranges 
from 20 to 50%, depending on pt. In the V2 vs mass method, the fraction 
of signal and background as well as the mass bin width were varied, giving 
a systematic uncertainty of about 30-40% pt-wide. 

The QC{2} method is sensitive to the variation of the background fit func- 
tion and to the D° selection cuts which give a total systematic uncertainty 
from 40% to 100% from low to high pt. 
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Fig. 3. Left panel: «2 vs pt with the 2A<fi bins method f ^2.2.1j) . boxes are systematic 
uncertainties. Right panel: V2 in 2 < p t < 5 GeV/c with all methods described. 



3. Results 

In Fig. [3j left panel, the D° V2 measured with the event plane method 
of the two A0 bins (full circle) for three pt intervals and the V2 for charged 
tracks measured by ALICE (empty star) are shown. The right panel shows 
the V2 from all methods in the first p t bin, where errors are the smallest, 
and they are consistent. The measurement is performed for inclusive D°, 
no feed-down subtraction was considered so far. 

Within 1.8cr a non-zero V2 is measured for D° mesons in the pt range 2-3 
GeV/c for Pb-Pb collisions at 2.6 TeV. 
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